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NATIONAL FOREWORD 

This Indian Standard (Part 2) which is identical with ISO 12213-2: 1997 'Natural gas— Calculation of 
compression factor — Part 2 : Calculation using molar-composition analysis' issued by the International 
Organization for Standardization (ISO) was adopted by the Bureau of Indian Standards on the 
recommendation of the Natural Gas Sectional Committee and approval of the Petroleum, Coal and Related 
Products Division Council. 

The text of International Standard has been approved as suitable for publication as Indian Standard 
without deviations. Certain conventions are, however, not identical to those used in Indian Standards. 
Attention is particularly drawn to the following: 

a) Wherever the words 'International Standard' appear referring to this standard, they should 
be read as 'Indian Standard'. 

b) Comma (,) has been used as a decimal marker while in Indian Standards, the current practice 
is to use a point (.) as the decimal marker. 

In this adopted standard reference appears to the following International Standards for some of which 
Indian Standards also exist. The corresponding Indian Standards, which are to be substituted in their 
place, are listed below along with their degree of equivalence for the editions indicated: 



International Standard 



Corresponding Indian Standard 



Degree of 
Equivalence 



ISO 31-3:1992 Quantities and 
units — Part 3 : Mechanics 

ISO 31-4:1992 Quantities and 
units — Part 4 : Heat 

ISO 6976:1995 Natural gas — 
Calculation of calorific values, 
density, relative density and 
Wobbe index from composition 

18012213-1:1997 Natural gas- 
Calculation of compression 
factor — Part 1 : Introduction and 
guidelines 



Nil 



Nil 



IS 14504:1998/180 6976:1995 Natural 
gas — Calculation of calorific values, 
density, relative density and Wobbe index 
from composition 

IS 15305 (Part 1):2003/ISO 12213-1:1997 
Natural gas— Calculation of compression 
factor:Part 1 Introduction and guidelines 



Identical 



do 



The Technical Committee responsible for the preparation of this standard has reviewed the provision of 
the ISO 31-3:1992 and ISO 31-4:1992 and decided that it is acceptable for use in conjunction with this 
standard. 
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Indian Standard 

NATURAL GAS — CALCULATION OF 

COMPRESSION FACTOR 

PART 2 CALCULATION USING MOLAR-COMPOSITION ANALYSIS 



1 Scope 

This International Standard specifies methods for the calculation of compression factors of natural gases, natural 
gases containing a synthetic admixture and similar mixtures at conditions under which the mixture can exist only as 
a gas. 

This part of ISO 12213 specifies a method for the calculation of compression factors when the detailed composition 
of the gas by mole fractions is known, together with the relevant pressures and temperatures. 

The method is applicable to pipeline quality gases within the ranges of pressure p and temperature T at which 
transmission and distribution operations normally take place, with an uncertainty of about ± 0,1 %. It can be applied, 
with greater uncertainty, to wider ranges of gas composition, pressure and temperature (see annex E). 

More detail concerning the scope and field of application of the method is given in part 1 of this International 
Standard. 



2 Normative references 

The following standards contain provisions which, through reference in this text, constitute provisions of this part of 
ISO 12213. At the time of publication, the editions indicated were valid. All standards are subject to revision, and 
parties to agreements based on this part of ISO 12213 are encouraged to investigate the possibility of applying the 
most recent editions of the standards indicated below. Members of lEC and ISO maintain registers of currently valid 
International Standards. 

ISO 31 -3: 1 992, Quantities and units — Part 3: Mechanics. 

ISO 31 -4: 1 992, Quantities and units — Part 4: Heat. 

ISO 6976:1995, Natural gas — Calculation of calorific values, density, relative density and Wobbe index from 
composition. 

ISO 1 221 3-1:1 997, Natural gas — Calculation of compression factor — Part 1: Introduction and guidelines. 



3 Definitions 

All definitions relevant to the use of this part of ISO 12213 are given in part 1. 
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4 Method of calculation 

4.1 Principle 

The method recommended uses an equation based on the concept that pipeline quality natural gas may be 
uniquely characterized for calculation of its volumetric properties by component analysis. This analysis, together 
with the pressure and temperature, are used as input data for the method. 

The method uses a detailed molar-composition analysis in which all constituents present in amounts exceeding a 
mole fraction of 0,000 05 should be represented. Typically, this includes all alkane hydrocarbons up to about Cy or 
Cg together with nitrogen, carbon dioxide and helium. 

For other gases, additional components such as water vapour, hydrogen sulfide and ethylene need to be taken into 
consideration (see reference [1] in annex G). 

For manufactured gases, hydrogen and carbon monoxide are also likely to be significant components. 

4.2 The AGA8-92DC equation 

The compression factor is determined using the AGA8 detailed characterization equation (denoted hereafter as the 
AGA8-92DC equation). This is an extended virial-type equation. The equation is described in AGA Report No. Sli). it 

may be written as 



Z - 1 + fip^ - A f,c: + £c;(fo, - c^k^p!"" ) p?" exp(-c„p*'' ) ... (1) 



18 58 

:13 n = 13 

where 



Z is the compression factor; 

a is the second virial coefficient; 

Pm is the molar density (moles per unit volume); 

Pr is the reduced density; 

h,„ c,„ k„ are constants (see table B.1); 

C* are coefficients which are functions of temperature and composition. 

The reduced density pr is related to the molar density pm by the equation 

Pr - ^^Ar, • ■ ■ (2) 

where Af is a mixture size parameter. 
The molar density can be written as 

Pm = pK^RT) ... (3) 

where 

p is the absolute pressure; 
R is the universal gas constant; 
7 is the absolute temperature. 

Z is calculated as follows: first the values of B and C* (« = 13 to 58) are calculated, using relationships given in 
annex B. Equations (1) and (3) are then solved simultaneously for pm and Z by a suitable numerical method (see 
figure B.1). 
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4.3 input variables 



The input variables required for use with the AGA8-92DC equation are the absolute pressure, absolute temperature 
and molar connposition. 

The composition is required, by mole fraction, of the following components: nitrogen, carbon dioxide, argon, 
methane, ethane, propane, n-butane, methyl-2-propane (iso-butane), «-pentane, methyl-2-butane (iso-pentane), 
hexanes, heptanes, octanes, nonanes, decanes, hydrogen, carbon monoxide, hydrogen sulfide, helium, oxygen 

and water. 

NOTE — If the mole fractions of the heptanes, octanes, nonanes and decanes are unknown, then use of a composite Cq+ 
fraction may be acceptable. The user should carry out a sensitivity analysis in order to test whether a particular approximation 
of this type degrades the result. 

All components with mole fractions greater than 0,000 05 shall be accounted for. Trace components (such as 
ethylene) shall be treated as given in table 1 . 

If the composition is known by volume fractions, these shall be converted to mole fractions using the method given 
in ISO 6976. The sum of all mole fractions shall be unity to within 0,000 1 . 

4.4 Ranges of application 
4.4.1 Pipeline quality gas 

The ranges of application for pipeline quality gas are as defined below: 

absolute pressure f\/IPa « p =s12MPa 

temperature 263 K ^ r ^ 338 K 

superior calorific value 30 MJ m3 « Hq ^ 45 I^Jm3 

relative density 0,55 =£ ^ =£ 0,80 

The mole fractions of the natural-gas components shall be within the following ranges: 

methane 

nitrogen 

carbon dioxide 

ethane 

propane 

butanes 

pentanes 

hexanes 

heptanes 

octanes plus 

higher hydrocarbons 
hydrogen 
carbon monoxide 
helium 
water 

Any component for which x, is less than 0,000 05 can be neglected. 

Minor and trace components are listed in table 1. 



0,7 


'=s.xcH4 «"1.00 





^ XN2 « 0,20 





« j:C02 ^ 0.20 





^.tC2H6 ^0.10 





=Sxc3H8 « 0,035 





^^C4Hio^ 0,015 





^ AC5H12 ^ 0-005 





=s xce =£ 0,001 





== xcy =s 0,000 5 





^ ^C8+ ^ 0,000 5 


0^ 


SXH2 *^0,10 


■ 


e-vco ^0,03 





sSxHe ^0,005 





^^H20^ 0,000 15 
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Table 1 — Minor and trace components 



Minor or trace component 


Assigned component 


Oxygen 


oxygen 


Argon 


argon 


Hydrogen sulfide 


hydrogen sulfide 


Ethylene, acetylene 


carbon dioxide 


Propylene, propadiene 


propane 


Butanes, butadienes 


n-butane 


Neo-pentane, pentenes, benzene, cyclopentane 


n-pentane 


All Ce-isomers, cyclohexane, ethylbenzene, xylenes 


rt-hexane 


All Cy-isomers, cycloheptane, toluene 


«-heptane 


All Cg-isomers 


n-octane 


All Cg-isomers 


n-nonane 


All CiQ-isomers and all higher hydrocarbons 


n-decane 



OMPa 


=£p 


^ 65 MPa 


225 K 


« T 


^350K 


0,55 


^d 


^0,90 



The metliod applies only to mixtures in the single-phase gaseous state (above the dew point) at the conditions of 
temperature and pressure of interest. 

4.4.2 Wider ranges of application 

The ranges of application tested beyond the limits given in 4.4.1 are: 

absolute pressure 

temperature 

relative density 

superior calorific value 20 MJ m-3 ^ //s =s 48 MJ m-3 

The allowable mole fractions of the major natural gas components are: 

methane 

nitrogen 

carbon dioxide 

ethane 

propane 

hydrogen 

The limits for minor and trace gas components are as given in 4.4.1 for pipeline quality gas. For use of the method 
outside these ranges, see annex E. 



0,50 


«-tCH4 «^-1.00 





s A-N2 s; 0,50 





«XC02 «^0,30 





« -«C2H6 «^ 0,20 





« -tCaHs =s 0,05 





=£j:H2 =£0,10 



4,5 Uncertainty 

4.5.1 Uncertainty for pipeline quality gas 

The uncertainty of results for use on all pipeline quality gas within the limits described in 4.4.1 is ± 0,1 % (for the 
temperature range 263 K to 350 K and pressures up to 12 MPa) (see figure 1). For temperatures above 290 K and 
at pressures up to 30 MPa the uncertainty of the result is also ± 0,1 %. 



For lower temperatures, the uncertainty of ± 0,1 % is at least maintained for pressures up to about 10 MPa. 
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A6A6-DC92 equation 



A AZ < ±0,1 % 

B ±0,1 % to ±0,2 % 

C ±0,2% to ±0,5% 



250 



275 



300 



325 



350 
Temperature, T (K) 



Figure 1 — Uncertainty limits for the calculation of compression factors (The uncertainty limits given are 

expected to be valid for natural gases and similar gases with j:n2 ^ ^,20, xco2 ^ 0,20, j:c2H6 ^0,10 

and JCH2 s; 0,10, and for 30 MJ-m-3 « f/g ^ 45 MJm-3 and 0,55 « rf « 0,80) 



This uncertainty level has been determined by comparison with the GERG databank of measurements of the 
compression factor for natural gases [2], [3]. a detailed comparison was also made with the GRI PVTdata on 
gravimetrically prepared simulated natural-gas mixtures W. [5]. 

The uncertainty of the measurements in both databanks used to test the method is of the order of ± 0,1 %. 

4.5.2 Uncertainty for wider ranges of application 

The estimated uncertainties for calculations of compression factors beyond the limits of quality given in 4.4.1 are 
discussed in annex E. 



4.5.3 Impact of uncertainties of input variables 

Listed in table 2 are typical values for the uncertainties of the relevant input variables. These values may be 
achieved under optimum operating conditions. 

As a general guideline only, an error propagation analysis using the uncertainties in the input variables produces an 
additional uncertainty of about ±0,1 % in the result at 6 MPa and within the temperature range 263 K to 338 K. 
Above 6 MPa, the additional uncertainties are greater and increase roughly in direct proportion to the pressure. 
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Table 2 — Uncertainties of input variables 



Input variable 


Absolute uncertainty 


Absolute pressure 


± 0,02 MPa 


Temperature 


+ 0,15 K 


Mole fraction of 




inerts 


±0,001 


nitrogen 


±0,001 


carbon dioxide 


+ 0,001 


methane 


±0,001 


ethane 


±0,001 


propane 


+ 0,000 5 


butanes 


± 0,000 3 


pentanes plus higher hydrocarbons 


± 0,000 1 


hydrogen and carbon monoxide 


± 0,001 



4.5.4 Reporting of results 

Results for compression factor and molar density shall be reported to four and to five places of decimals, 
respectively, together with the pressure and temperature values and the calculation method used (ISO 12213-2, 
AGA8-92DC equation). For verification of calculation procedures, it is useful to carry extra digits. 



5 Suppliers of computer programmes 

It is planned to make software available which implements this International Standard. Users are invited to contact 
their ISO member body or ISO Central Secretariat to enquire about the availability of such software. 
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Annex A 

(normative) 

Symbols and units 



Symbol 



B 
B 



"U 



•- n 



F 

fn 

G 

C, 



G 



r . . 

'./ 

En 

Us 

K 
K, 

'<J 
'<<J 

kn 

M 

M, 
N 
n 
P 

Q 

a 

Qj 

R 



Meaning 

Constant in table B.I 
Second virial coefficient 

Mixture interaction coefficient [equations (B.1) and (B.2)] 

Constant in table B.I 
Constant in table B.1 

Coefficients whicii are functions of temperature and composition 

Characteristic energy parameter for (th component (table' B.2) 
Characteristic energy parameter for yth component (table B.2) 
Binary energy parameter for second virial coefficient 

Binary energy interaction parameter for second virial coefficient (table B.3) 

Mixture high-temperature parameter 

High-temperature parameter for /th component (table B.2) 

High-temperature parameter f or ;th component (table B.2) 

Constant in table B.1 

Mixture orientation parameter 

Orientation parameter for Hh component (table B.2) 

Orientation parameter for jth component (table B.2) 

Binary orientation parameter 

Binary interaction parameter for orientation (table B.2) 

Constant in table B.1 

Superior calorific value 

Size parameter 

Size parameter for /th component (table B.2) 

Size parameter foryth component (table B.2) 

Binary interaction parameter for size (table B.3) 

Constant in table B.1 

Molar mass 

Molar mass of I'th component 

Number of components in gas mixture 

An integer (from 1 to 58) 

Absolute pressure 

Quadrupole parameter 

Quadrupole parameter for ith component 

Quadrupole parameter foryth component 

Constant (table B.1) 

Gas constant (= 0,008 314 510) 

Dipole parameter for ith component (table B.2) 

Dipole parameter foryth component (table B.2) 

Constant (table B.1) 



Units 

m3 kmol-i 



K 
K 
K 



MJm3 
(m3/kmol)i/3 
(m3/kmol)i/3 
(m3/kmol)i/3 



kgkmol-i 
kgkmoH 



MPa 



MJ(kmolK)- 
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Symbol Meaning Units 

T Absolute temperature K 

U Mixture energy parameter K 

Uij Binary interaction parameter for mixture energy (table B.3) — 

M„ Constant in table B.I — 

W, Association parameter for ith component (table B.2) — 

Wj Association parameter foryth component (table B.2) — 

M„ Constant (table B.1) — 

Xj Mole fraction of iVc\ component in gas mixture — 

Xj Mole fraction of yth component in gas mixture — 

Z Compression factor — 

p Mass density kgm-3 

Pr Reduced density of gas — 

Pm Molar density l<molm-3 
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Annex B 

(normative) 

Description of the AGA8-92DC method 



B.1 General 



For gas mixtures, the compression factor Z is calculated using the equations given in 4.2. This annex gives a 
detailed description ot the computations and the necessary numerical values. The description is based upon that 
given in AGA Report No. 8 [1]. A programme implementing this description is given in annex F, and as such 
provides the correct solution. Other computational procedures are acceptable provided that they can be 
demonstrated to yield identical numerical results {see annex C for examples). 



B.2 Computer implementation of the AGA8-92DC method 
B.2.1 Overview of the calculation procedure 

I Input the absolute temperature T, absolute pressure p and mole fraction of each component x; of the mixture. 

NOTE — For pressure and temperature, values known in any other units will first have to be converted precisely to 
values In megapascals and kelvlns, respectively (see ISO 31-3 and ISO 31-4 and annex D for relevant conversion 
factors). 

II Compute the equation of state coefficients B and C*„ (n = 13 to 58) that depend on Tand Xi. 

III Solve iteratively for the molar density p^, using the equation of state rearranged to give the pressure p. 

IV Output the compression factor after the computed pressure from step III and the input pressure from step I 
agree within a specified range of convergence (e.g. 1E-06). 

Figure B.I shows a flow diagram of these steps. 

B.2.2 Details of the calculation procedure 
Step! 

Input the absolute temperature T, absolute pressure p and mole fraction a-, of each constituent in the natural-gas 

mixture. 

Step II 

At the absolute temperature T and the mole fractions xy of the natural gas (as input from step I), compute the 
composition- and temperature-dependent coefficients B and C^ (n = 13 to 58). 

The second virial coefficient B is given by the following equations: 

18 N N 3-2 

« = 1 / = 1 y = 1 

fi;,^^(G,^ +i_g„)^'"(e.e. +^^q^y''(F^^FJ^^^-f„f(s^sJ^^-s,)'"{w^wJ -n-w„p . ..(b.2) 
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Input T. p. X/ 



Compute 

equation of state coefficients 

B and C'„ (n = 13 to 58) 




Figure B.I — AGA8-92DC equation — Calculation flow diagram 



The binary parameters Eij and Cy are calculated using the following equations: 

Efj = E;j(E,Ejf^ ...(B.3) 

G,j=G;j(Gi+Gj)/2 ...(B.4) 

Note that all values of the binary interaction parameters £"* and G* are 1 ,0 except for the values given in table B.3. 
The coefficients C* (n = 13 to 58) are given by the equation: 

C; =a„{G + ^- g, Y" (^2 + 1 - <7„ )''" {F + ^-f„ Y" U"n T-"n . . . (B.5) 

The mixture parameters V, G and Q are calculated using the following conformal solution mixing equations, where 
in the double sums / ranges from 1 to A? - 1 and, for each value of i, j ranges from i + 1 to A': 



i/5 = 



X:E, 



5/2 



W-1 yv 
1=1 ;=i+1 



5/2 



N N -■[ f^ 

(=1 ,=1 y=,+1 



. . . (B.6) 



(B.7) 



10 



1=1 

N 
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. . . (B.8) 



(8.9) 



It should be noted that all values of the binary interaction parameters A:,y, E'^j , G* and f/y are 1,0 except for the 
values given in table B.3. Also note that F, is zero for all components except hydrogen, for which F(H2) = 1,0, and 
that W, is zero for all components except water, for which IV(H20) = 1 ,0. 



Step III 

In the computation of the compression factor Z, the composition of the gas is known, the absolute temperature T of 
the gas is known and the absolute pressure is known. The problem then is to compute the molar density pm, using 
the equation of state expression for the pressure p. For this purpose, the definition of the compression factor Z as 
given in equation (1) (see 4.2) is substituted into equation (3) to obtain an equation for the pressure as given in 
equation (B.10): 



P = PmRT 



18 58 , ic \ h I k 



n = ^3 



/i = 13 



(B.10) 



Equation (B.10) is solved using standard equation of state density search algorithms. Having obtained an equation 
for the pressure p [equation (B.10)], the problem is then to search for the value of the molar density pm that will yield 
the pressure that is within a preset limit (e.g. 1 x 10-6) equal to the input pressure. 

The reduced density pr is related to the molar density pm by the mixture size parameter [see equation (2) in 4.2], 

The mixture size parameter K is calculated using the following equation: 



K5 = 



f. 



N-^ 



1=1 7 =1+1 



5/2 



(B.11) 



Note that in the summations the subscript / refers to the /th component in the gas mixture and the subscript 7 refers 
to theyth component in the gas mixture. The quantity N is the number of components in the mixture. Thus, in the 
single summation, i ranges over the integer values from 1 to N. For example, for a mixture of 12 components, 
N= 12 and there would be 12 terms in the single sum. In the double summation, / ranges from 1 to A'- 1 and, for 
each value of /, j ranges from j + 1 to N. For example, for a mixture of 12 components, there would be 66 terms in 
the double summation if all values of Kij differed from 1,0. However, because many of the values of Kij are 1,0, the 
number of non-zero terms in the double summation is small for many natural-gas mixtures. Note that all values of 
Kij are 1 ,0 except for the values given in table B.3. 

Step IV 

Once the molar density p^ has been obtained in step III, the compression factor is calculated in step IV using the 
pressure, temperature, molar density and gas constant: 



NOTE — The density p (mass per unit volume) can be calculated as follows: 

P = MPtn 



.(8.12) 



(B.13) 



11 
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where M is calculated from the equation: 



M 



(B.14) 



Report the density to three places of decimals. 

Table B.I — Equation of state parameters 



n 


"n 


bn 


c„ 


kn 


"n 


Sn 


<?„ 


fn 


Sn 


»'n 


1 


0,153 832 600 










0,0 

















2 


1 ,341 953 000 










0.5 

















3 


- 2,998 583 000 










1,0 

















4 


-0,048 312 280 










3.5 

















5 


0,375 796 500 










-0,5 


1 














6 


- 1 ,589 575 000 










4,5 


1 














7 


- 0,053 588 470 










0.5 





1 











8 


0,886 594 630 










7,5 











1 





9 


-0.710 237 040 










9,5 











1 





10 


-1,471 722 000 










6,0 














1 


11 


1,321 850 350 










12,0 














1 


12 


- 0,786 659 250 










12,5 














1 


13 


2,291 290 X 10-9 






3 


-6,0 








1 








14 


0,157 672 400 






2 


2,0 

















15 


- 0,436 386 400 






2 


3,0 

















16 


- 0,044 081 590 






2 


2,0 





1 











17 


- 0,003 433 888 






4 


2,0 

















18 


0,032 059 050 






4 


11,0 

















19 


0,024 873 550 


2 







-0,5 

















20 


0,073 322 790 


2 







0,5 

















21 


- 0,001 600 573 


2 




2 


0,0 

















22 


0,642 470 600 , 


2 




2 


4,0 

















23 


-0,416 260 100 


2 




2 


6,0 

















24 


- 0,066 899 570 


2 




4 


21,0 

















25 


0,279 179 500 


2 




4 


23,0 


1 














26 


-0,696 605 100 


2 




4 


22,0 





1 











27 


- 0,002 860 589 


2 




4 


-1,0 








1 








28 


- 0,008 098 836 


3 







-0,5 





1 











29 


3,150 547 000 


3 




1 


7,0 


1 














30 


0,007 224 479 


3 




1 


-1,0 








1 








31 


- 0,705 752 900 


3 




2 


6,0 

















32 


0,534 979 200 


3 




2 


4.0 


1 














33 


-0,079 314 910 


3 




3 


1,0 


1 














34 


-1,418 465 000 


3 




3 


9,0 


1 














35 


-5,999 05x10-17 


3 




4 


-13,0 








1 








36 


0,105 840 200 


3 




4 


21,0 

















37 


0,034 317 290 


3 




4 


8,0 





1 











38 


- 0,007 022 847 


4 







-0,5 

















39 


0,024 955 870 


4 







0,0 

















40 


0,042 968 180 


4 




2 


2,0 

















41 


0,746 545 300 


4 




2 


7,0 

















42 


- 0,291 961 300 


4 




2 


9,0 





1 











43 


7,294 616 000 


4 




4 


22,0 

















44 


- 9,936 757 000 


4 




4 


23,0 

















45 


- 0,005 399 808 


5 








1,0 

















46 


- 0,243 256 700 


5 




2 


9,0 

















47 


0,049 870 160 


5 




2 


3,0 





1 











48 


0,003 733 797 


5 




4 


8,0 

















49 


1 ,874 951 000 


5 




4 


23,0 





1 











50 


0,002 168 144 


6 








1,5 

















51 


-0,658 716 400 


6 




2 


5,0 


1 














52 


0,000 205 518 


7 








-0,5 





1 











53 


0,009 776 195 


7 




2 


4,0 

















54 


- 0,020 487 080 


8 




1 


7,0 


1 














55 


0,015 573 220 


8 




2 


3,0 

















56 


0,006 862 415 


8 




2 


0,0 


1 














57 


- 0,001 226 752 


9 




2 


1,0 

















58 


0,002 850 908 


9 




2 


0,0 





1 
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Table B.2 — Characterization parameters 



identifi- 
cation 
number 



Compound 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 



Methane 

Nitrogen 

Carbon dioxide 

Ethane 

Propane 

Water 

Hydrogen sulfide 

Hydrogen 

Carbon monoxide 

Oxygen 

iso-Butane 

n- Butane 

iso-Pentane 

n-Pentane 

n-Hexane 

n-Heptane 

«-Octane 

n-Nonane 

n-Decane 

Helium 

Argon 



Molar mass 



kgkmol-i 

16,043 

28,013 5 

44,010 

30,070 

44,097 

18,0153 

34,082 

2,015 9 

28,010 

31,998 8 

58,123 

58,123 

72,150 

72,150 

86,177 

100,204 

114,231 

128,258 

142,285 

4,002 6 

39,948 



Energy 
parameter 



K 

151,318 300 
99,737 780 
241,960 600 
244,166 700 
298,118 300 
514,015 600 
296,355 000 
26,957 940 
105.534 800 
122,766 700 
324.068 900 
337,638 900 
365,599 900 
370,682 300 
402,636 293 
427,722 630 
450,325 022 
470,840 891 
489,558 373 
2,610 111 
119,629 900 



Size parameter Orientation Quadrupole High- Dipole Associ- 
parameter parameter temp, parameter ation 
param- param- 
eter eter 



(m3/kmol)'"3 

0,461 925 5 
0,447 915 3 
0,455 748 9 
0,527 920 9 
0,583 749 
0,382 586 8 
0,461 826 3 
0,351 491 6 
0,453 389 4 
0,418 695 4 
0,640 693 7 
0,634 142 3 
0,673 857 7 
0,679 830 7 
0,717 511 8 
0,752 518 9 
0,784 955 
0,815 273 1 
0,843 782 6 
0,358 988 8 
0,421 655 1 



0,0 

0,027 815 
0,189 065 
0,079 300 
0,141 239 
0,332 500 
0,088 500 
0,034 369 
0,038 953 
0,021 000 
0,256 692 
0,281 835 
0,332 267 
0,366 911 
0,289 731 
0,337 542 
0,383 381 
0,427 354 
0,469 659 
0,0 
0,0 



Q. 



0,0 

0,0 

0,690 000 

0,0 

0,0 

1 ,067 750 

0,633 276 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0,0 

0.0 

0,0 

0,0 

0,0 

0,0 



0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 1 .582 200 

0,0 0,390 000 

1,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 

0,0 0,0 



W: 



0,0 
0,0 
0,0 
0,0 
0,0 
1,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 
0,0 



Table B.3 — Binary interaction parameter values 



Identification 


Component pair 


Eij 


^ij 


l^ij 


C,> 


nunnber 














'■ J 
1 2 


Methane - 


1- nitrogen 


0.971 640 


0,886 106 


1,003 630 




3 




carbon dioxide 


0,960 644 


0,963 827 


0,995 933 


0,807 653 


4 




ethane 










5 




propane 


0,994 635 


0,990 877 


1,007 619 




6 




water 


0,708 218 








7 




hydrogen sulfide 


0,931 484 


0,736 833 


1 ,000 080 




8 




hydrogen 


1,170 520 


1,156 390 


1 ,023 260 


1,957 310 


9 




carbon monoxide 


0,990 126 








10 




oxygen 










11 




iso-butane 


1,019 530 








12 




n-butane 


0,989 844 


0,992 291 


0,997 596 




13 




iso-pentane 


1,002 350 








14 




fi-pentane 


0,999 268 


1 ,003 670 


1 ,002 529 




15 




n-hexane 


1,107 274 


1,302 576 


0,982 962 




16 




«-heptane 


0,880 880 


1,191 904 


0,983 565 




17 




n-octane 


0,880 973 


1 ,205 769 


0,982 707 




18 




n-nonane 


0,881 067 


1,219 634 


0,981 849 




19 




«-decane 


0,881 161 


1,233 498 


0,980 991 
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Identification 


Component pair 


E, 


^'J 


f^>j 


G, 


number 














i 
2 


J 
3 


Nitrogen + 


carbon dioxide 


1,022 740 


0,835 058 


0,982 361 


0,982 746 




4 




ethane 


0,970 120 


0,816 431 


1 ,007 960 






5 




propane 


0,945 939 


0.915 502 








6 




water 


0,746 954 










7 




hydrogen sulfide 


0,902 271 


0,993 476 


0,942 596 






8 




hydrogen 


1,086 320 


0,408 838 


1,032 270 






9 




carbon monoxide 


1,005 710 










10 




oxygen 


1,021 000 










11 




iso-butane 


0,946 914 










12 




n-butane 


0,973 384 


0,993 556 








13 




iso-pentane 


0,959 340 










14 




n-pentane 


0,945 520 








3 


4 
5 


Carbon dioxide + 


ethane 
propane 


0,925 053 
0,960 237 


0,969 870 


1,008 510 


0,370 296 




6 




water 


0,849 408 






1,673 090 




7 




hydrogen sulfide 


0,955 052 


1.045 290 


1,007 790 






8 




hydrogen 


1 ,281 790 










9 




carbon monoxide 


1,500 000 


0,900 000 








10 




oxygen 












11 




iso-butane 


0,906 849 










12 




n-butane 


0,897 362 










13 




iso-pentane 


0,726 255 










14 




,1-pentane 


0,859 764 










15 




n-hexane 


0,855 134 


1 ,066 638 


0,910 183 






16 




ra-heptane 


0,831 229 


1,077 634 


0,895 362 






17 




n -octane 


0,808 310 


1,088 178 


0,881 152 






18 




n-nonane 


0,786 323 


1,098 291 


0,867 520 






19 




tt-decane 


0,765 171 


1,108 021 


0,854 406 




4 


5 
6 


Etiiane + 


propane 
water 


1 ,022 560 
0,693 168 


1,065 173 


0,986 893 






7 




hydrogen sulfide 


0,946 871 


0,971 926 


0,999 969 






8 




hydrogen 


1,164 460 


1,616 660 


1,020 340 






9 




carbon monoxide 












10 




oxygen 












11 




iso-butane 




1 ,250 000 








12 




n-butane 


1,013 060 


1 ,250 000 








13 




iso-pentane 




1 ,250 000 








14 




n-pentane 


1,005 320 


1 ,250 000 






5 


8 
12 


Propane + 


hydrogen 
n-butane 


1.034 787 
1 ,004 900 








7 


15 


Hydrogen sulfide + 


n-hexane 


1,008 692 


1,028 973 


0,968130 






16 




«-heptane 


1,010 126 


1 ,033 754 


0.962870 






17 




n-octane 


1,011 501 


1,038 338 


0,957828 






18 




Ai-nonane 


1,012 821 


1,042 735 


0,952 441 






19 




«-decane 


1,014 089 


1 ,046 966 


0.948 338 




8 


9 
10 
11 
12 


Hydrogen + 


carbon monoxide 
oxygen 
Iso-butane 
n-butane 


1,100 000 

1 ,300 000 
1 ,300 000 
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Annex C 

(normative) 

Example calculations 

The following example calculations have been carried out using the validated computer programme described in 
reference [1], which incorporates the subroutine described in annex B. 

Table C.I — Gas analysis in mole fractions 





Gasi 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


*C02 


0,006 


0,005 


0,015 


0,016 


0,076 


0,011 


XN2 


0,003 


0,031 


0,010 


0,100 


0,057 


0,117 


^H2 


0,00 


0,00 


0,00 


0,095 


0,00 


0,00 


^CO 


0,00 


0,00 


0,00 


0,010 


0,00 


0.00 


^CH4 


0,965 


0,907 


0,859 


0,735 


0,812 


0,826 


*C2H6 


0,018 


0,045 


0.085 


0,033 


0,043 


0,035 


•«C3H8 


0,004 5 


0,008 4 


0,023 


0.007 4 


0,009 


0,007 5 


'Ciso-C4HiO 


0,001 


0.001 


0,003 5 


0,001 2 


0,001 5 


0,001 2 


Xn-CAH-\0 


0,001 


0,001 5 


0,003 5 


0,001 2 


0,001 5 


0,001 2 


J^iso-C5Hl2 


0,000 5 


0,000 3 


0,000 5 


0,000 4 


0,00 


0,000 4 


->^n-C5Hl2 


0,000 3 


0,000 4 


0,000 5 


0,000 4 


0,00 


0,000 4 


•»^C6H14 


0,000 7 


0,000 4 


0,00 


0,000 2 


0,00 


0,000 2 


■«C7Hi6 


0,00 


0,00 


0,00 


0,000 1 


0,00 


0,000 1 


J^CsHls 


0,00 


0,00 


0,00 


0,000 1 


0,00 


0,00 



Table C.2 — Results (Z-values) 



Cond 


itions 


Gasi 


Gas 2 


Gas 3 


Gas 4 


Gas 5 


Gas 6 


P 
bar 


°c 














60 


-3,15 


0,840 53 


0,833 48 


0,793 80 


0,885 50 


0,826 09 


0,853 80 


60 


6,85 


0,861 99 


0,855 96 


0,822 06 


0,901 44 


0,849 69 


0,873 70 


60 


16,85 


0,880 06 


0,874 84 


0,845 44 


0,915 01 


0,869 44 


0,890 52 


60 


36,85 


0,908 67 


0,904 46 


0,881 83 


0,936 74 


0,900 52 


0,917 23 


60 


56,85 


0,930 1 1 


0,926 96 


0,908 68 


0,953 18 


0,923 68 


0,937 30 


120 


-3,15 


0,721 33 


0,710 44 


0,641 45 


0,810 24 


0,695 40 


0,750 74 


120 


6,85 


0,760 25 


0,750 66 


0,689 71 


0,837 82 


0,737 80 


0,785 86 


120 


16,85 


0,793 17 


0,784 75 


0,731 23 


0,861 37 


0,773 69 


0,815 69 


120 


36,85 


0,845 15 


0,838 63 


0,796 98 


0,899 13 


0,830 22 


0,863 1 1 


120 


56,85 


0,883 83 


0,878 70 


0,845 53 


0.927 66 


0,872 1 1 


0,898 62 
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Annex D 

(normative) 

Pressure and temperature conversion factors 

If the input variables pressure and temperature are not in the necessary units of megapascals and Kelvins, then 
conversions will have to be made in order to use the Fortran implementation. A selection of appropriate conversion 
factors is given below: 

Pressure 

/5(MPa) = [p(bar)]x10-i 

p (MPa) = b(atm)] x 0,101 325 

/p(MPa) = [;;(psia)]/1 45,038 

P (MPa) = b(psig) + 14,695 9]/1 45,038 

Temperature 

r(K) = f{°C) + 273,15 

T (K) = [/(°F) - 32]/1 ,8 + 273,15 

7(K) = [rrR)]/1,8 
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Annex E 

(informative) 

Performance over wider ranges of application 

The AGA8-92DC equation has been comprehensively tested, over the temperature range 263 K to 338 K and at 
pressures up to 30 MPa, with the GERG data [2]. [31 and the Gas Research Institute data W for gases within the 
ranges of composition given for pipeline quality gases (see 4.4.1). Within these limits, the uncertainties are as given 
in 4.5. 

Rough estimates of the uncertainties involved in calculations of compression factors for wider ranges of application 
(with respect to composition) (see 4.4,2) are plotted in figures E.I to E.4 as pressure-composition plots for nitrogen, 
carbon dioxide, ethane and propane, respectively. 

In figures E.I to E.4 the performance of the AGA8-92DC method is illustrated up to a maximum pressure of 
30 MPa. The uncertainty limits are dependent upon pressure, temperature and composition, and are also strongly 
affected by the proximity of the phase boundary. The estimated uncertainty limits presented below are based upon 
less comprehensive data, published as a supplement to the GERG databank PI, and upon the databanks in 
references [2] and [4]. Reference [4] also provides data up to 70 MPa. The uncertainty limits given in figures E.1 to 
E.4 are always for the worst-case result, i.e. they are the least optimistic choice. 

Dashed lines are used to separate two regions of estimated uncertainty when the experimental evidence is not 
sufficient to determine the position of the boundary. The detailed composition of the gas will have a strong influence 
on the position of the phase boundary and the user should, therefore, make his own phase boundary calculation. 




AGAe-0C92 equation 

r= 263 K to 338 K 



0,2 



0,3 



0,4 0,5 

Mole fraction of nitrogen, x^, 



Figure E.I — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high nitrogen content 
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- 30 



20 



10 



AGAe-DC92 equation 

r = 263K to 338 K 



A AZ 5 iO,1 % 

B ±0,1% to ±0,2% 

C ±0,2% to ±0,5% 

D ±0,5% to ±3,0% 




0,1 



0,2 0,3 

Mole fraction of carbon dioxide, x^a? 



Figure E.2 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high carbon dioxide content 
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30 



20 



10 




h AZ < ±0,1 % 

B ±0,1 % to ±0,2 % 

C ±0,2% to ±0,5% 



0,15 0,2 

Mole fraction of etfiane, x^2hc 



Figure E.3 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high ethane content 
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2; 



20 



10 



AGA8-DC92 equation 

7" = 263 K to 338 K 



A A/ £ ±0,1 % 

B ±0,1% to ±0,2% 



0,05 



0,1 
Mole fraction of propane, >•[ 



Figure E.4 — Estimated uncertainty limits for the calculation of compression factors of natural gases with a 

high propane content 

The overall results at pressures up to 10MPa and temperatures within the range 263 K to 338 K can be 
sunnmarized as follows. Only gases having mole fractions within the limits given below will have uncertainties within 
± 0,1 %, ± 0,2 % and ± 0,5 %, respectively, within the given pressure and temperature domain. 



Component 


Mole fraction for an uncertainty within | 




± 0,1 % 


+ 0,2 % 


± 0,5 % 


Nitrogen 


s;0,50 


— 


— 


Carbon dioxide 


=s0,23 


=5 0,26 


^0,28 


Ethane 


«0,13 


«0,20 


— 


Propane 


=s0,06 


«0,10 


— 
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Annex F 

(informative) 

Subroutines in Fortran for thie AGA8-92DC method 



c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
r 



update: 
update : 



17.05.94 
4.09.95 



E.W. Lemmon/S.W. Beyerlein/ J . L . Savidge 
M. Jaeschke J. Sikora 



AGA8-DC92 COMPRESSION FACTOR EQUATION 

SUBROUTINE DCAGA 

This program was written to accompany ISO 12213. 

"DCAGA" Calculates the compression factor of natural gases using 
a detailed gas analysis. 

For information contact: DR. Jeffrey L. Savidge 

gas research institute 
8 600 W. Bryn Mawr Ave. 
Chicago, IL 60631 
(312) 399-8100, FAX (312) 399-8125 

This program calculates compression factors and molar densities for 
natural gases from the input of gas composition in accordance with the 
AGA8-DC92 compression equation developed by the Gas 
Research Institute, Chicago, Illinois. (K.E. Starling and J.L. 
Savidge, Compressibility Factors of Natural Gas and Other Related 
Hydrocarbon Gases, American Gas Association, AGA Transmission 
Measurement Committee Report No. 8, American Petroleum Institute MPMS 
Chapter 14.2, Second Edition, November 1992, Catalog No. XQ9212). 

The coefficients used in this program are the same as the values found 
in AGA Report No. 8, November 1992. 

Values for the gas constant and molar masses conform with ISO 6976 
(1995) and GPA 2172 (1988) . 



Ranges of application for compression factor calculation 
with the AGA8 - DC92 equation: 



p-T-ranges 



absolute pressure 
temperature 



to 65 MPa 
225 to 350 K 



Ranges for percentage molar composition: 

A: pipeline quality natural gas 
B: wider ranges of application 



methane 

nitrogen 

carbon dioxide 

ethane 

propane 

butanes 

pentanes 

hexanes 

heptanes 

octanes plus 

hydrogen 

carbon monoxide 

helium 



70 to 100 



B 
50 to 100 



to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 



20 
20 
10 

3. 

1, 



0, 

0.05 

0.05 
10 

3 

0.5 



to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 
to 



50 

30 

20 
5 

1.5 
0.5 
0.1 
0.05 
0.05 

10 
3 
0.5 
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C water to 0.015 to 0.015 

C 

C The expected uncertainty of the calculated results are for 

C pipeline quality natural gases: 

C +/- 0,1% within p-T range to 12 MPa, 263 to 350 K 

C +/- 0,6% within p-T range to 60 MPa, 225 to 350 K 

C The expected uncertainty in the wider range of application 

C (composition) is often even for pressures up to 12 MPa larger. 

C For more details see ISO 12213 part 2 . 

C 

SUBROUTINE DCAGA (XJ) 

INTEGER B(58) ,C{58) ,K{58) ,G(58) 
INTEGER Q(58) ,F(58) ,S(58) ,W(58) 
REAL*8 A{58) ,U(58) 
COMMON /CONSTANTS/ A, B, C , K, U, G, Q, F, S, W 

REAL* 8 MW(21) ,EI{21) ,KI{21) ,GI(21) ,QI(21) ,FI{21) ,SI(21) ,WI(21) 

REAL* 8 EIJ{21,21) ,UIJ(21,21) ,KIJ(21,21) ,GIJ(21,21) 

COMMON /PARAMETERS/ MM, EI , KI , GI , QI , FI , SI , WI , EIJ, UIJ, KIJ, GIJ 

REAL*8 Kl, CNS{58), BI(18) 
COMMON /COEF/ Kl , CNS , BI 

REAL* 8 MWX, RGAS, TCM, DCM 
COMMON /MW/ MWX, RGAS, TCM, DCM 

INTEGER I, J, N 
REAL*8 SUM, XI (21), XJ(21) 
REAL* 8 Ul, Gl, Ql, Fl, El 
REAL*8 XIJ, EIJO, GIJO, BN 



XI { 1) = 


XJ( 1} 


XI ( 4) = 


XJ( 2) 


XI( 5) = 


XJ( 3) 


XKll) = 


XJ( 4) 


XI (12) = 


XJ{ 5) 


XI (13) = 


XJ( 6) 


XI(14) = 


XJ( 7) 


XI(15) = 


XJ( 8) 


XI(16) = 


XJ( 9) 


XI(17) = 


XJ(IO) 


XI(18) = 


XJ(ll) 


XI(19) = 


XJ<12) 


XI ( 3) = 


XJ(13) 


XI ( 2) = 


XJ(14) 


XI ( 7) = 


XJ(15) 


XI(20) = 


XJ(16) 


XI ( 6) = 


XJ(17) 


XI(IO) = 


XJ(18) 


XI(21) = 


XJ(19) 


XI( 8) = 


XJ(20) 


XI ( 9) = 


X J ( 2 1 ) 


Normalize molp 


SUM = 


.( 


MWX = 





DO 10 1=1, 21 
10 SUM = SUM + XI(I) 
DO 20 1=1, 21 

2 XI (I) = XI(I) /SUM 

C Calculate molecular weight 

RGAS = 8.31451D-3 

MWX = 

DO 30 1=1, 21 

3 MWX = MWX + XI(I)*MVJ(I) 

DO 4 N=l, 18 
40 BI(N) = 
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Kl = 

Ul = 

Gl = 

Ql = 

Fl = 

El = 

DO 50 1=1, 21 

Kl = Kl + XI(I) *KI(I) **2 .5D0 

Ul = Ul + XI(I) *EI(I) **2 .5D0 

Gl = Gl + XI(I)*GI(I) 

Ql = Ql + XI(I)*QI(I) 

Fl = Fl + XI(I)*XI(I) *FI(I) 

El = El + XI(I) *EI (I) 

5 CONTINUE 

TOM = 1.251*E1 
DCM = K1**(-1.2D0) 
Kl = K1*K1 
Ul = U1*U1 

DO 50 1=1, 8 

DO 60 J=I+1, 19 

XIJ = XI(I) *XI{J) 
IF (XIJ.NE.O) THEN 

Kl = K1+2.D0*XIJ*(KIJ(I, J)**5.D0-1.D0)* (KI(I)*KI(J)) **2.5D0 
Ul = U1+2.D0*XIJ*(UIJ{I, J) **5.DG-1.D0) *{EI(I) *EI(J) ) **2.5D0 
Gl = G1+XIJ*{GIJ(I,J) - l.DO)*(GI(I) + GI{J)) 
END IF 

6 CONTINUE 

DO 80 1=1, 21 
DO 80 J=I, 21 

XIJ = XI (I) *XI(J) 
IF (XIJ.NE.O) THEN 

IF (I.NE.J) XIJ = 2.D0*XIJ 
EIJO = EIJ(I, J) *DSQRT(EI (I) *EI(J) ) 
GIJO = GIJ(I, J)*(GI(I) + GI(J))/2.D0 
DO 70 N=l, 18 



'^F(N) 





BN 


= 


(GIJO 


+ l.DO 


- G 


(N) ) 


**G(N) 




So 




* 


(QKI) 


*QI(J) 


+ 1 


.DO 


- Q(N))**Q(N) 


& 




* 


(DSQRT(FI (I) * 


'FI(J) ) 


+ l.DO - 


F(N) ) 


& 




* 


(SKI) 


*SI(J) 


+ 1, 


.DO 


- S(N) ) **S(N) 


& 




* 


(WI(I) 


*WI(J) 


+ 1, 


.DO 


- W(N) ) **W(N) 




BI (N) 


= BI(N)+A(N) 


*XIJ*EIJO**U(N) ■ 


MKI(I 




CONTINUE 














ENDIF 
















CONTINUE 

















7 
80 

Kl = K1**0.2D0 

Ul = U1**0.2D0 

DO 90 N=13, 58 
90 CNS(N) = (Gl + l.DO - G(N))**G(N) 

& * (Ql**2 + l.DO - Q(N))**Q(N) 
& * (Fl + l.DO - F(N))**F(N) 
& * A(N) *U1**U(N) 

END 

SUBROUTINE PZOFDT(D, T, P, Z, BMIX) 

INTEGER B(58) ,C(58) ,K(58) ,G(58) 
INTEGER Q(58) ,F(58) ,S(58) ,W(58) 
REAL*8 A(58) ,U(58) 
COMMON /CONSTANTS/ A, B, C, K, U, G, Q, F, S , W 

REAL* 8 Kl, CNS(58), BI(18) 
COMMON /COEF/ Kl , CNS , BI 

REAL* 8 MWX, RGAS, TCM, DCM 
COMMON /MW/ MWX, RGAS, TCM, DCM 



5D0*BN 
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INTEGER N 

REAL* 8 D, T, P, Z, BMIX, DR 

DR = D*K1**3 
BMIX = 
DO 10 N=l, 18 
10 BMIX = BMIX + BI{N) /T**U(N) 

Z = l.DO + BMIX*D 
DO 20 N=13, 18 
20 Z = Z - DR*CNS(N) /T**U(N) 

DO 30 N=13, 58 
30 Z = Z + CNS(N) /T**U(N) * (B(N) - C (N) *K (N) *DR**K (N) ) *DR**B (N) 
& *DEXP(-C(N) *DR**K(N) ) 

P = D*RGAS*T-^2 
END 

SUBROUTINE DZOFPT(P, T, D, Z, BMIX) 

REAL*8 P, T, D, Z, BMIX 

REAL* 8 XI, X2, X3 , F, Fl, F2 , F3 , TOL 

INTEGER I 

TOL = 0.5D-9 
XI = O.OOOOOIDO 
X2 = 40.D0 
D = 



CALL PZ0FDT(X1, 


T, Fl, Z, 


BMIX) 


CALL PZ0FDT(X2, 


T, F2, Z, 


BMIX) 


Fl = Fl - P 






F2 = F2 - P 






IF (Fl*F2.GE.O) 


RETURN 





C BEGIN ITERATING 
C 



DO 60 I = 1, 50 
C ...Use False Position to get point 3. 
X3 = XI - F1*{X2 - XI) /(F2 - Fl) 
CALL PZ0FDT(X3, T, F3 , Z, BMIX) 
P3 = F3 - P 

C ...Use points 1, 2, and 3 to estimate the root using Chafhber ' s 
C . . .method (quadratic solution) . 
D = X1*F2*F3/ { (Fl 
& + X2*F1*F3/ ( (F2 
& + X3*F1*F2/ ( (F3 

IF ( (D - X1)*(D - X2).GE.O) D = (XI + X2)/2.D0 

CALL PZOFDT(D, T, 

F = F - P 

IF (DABS(F) .LE.TOL) 

C ...Discard quadratic solution if false position root is closer. 
IF (DABS(F3) .LT.DABS(F) .AND. F*F3.GT.0) THEN 
IF (F3*F1.GT.0) THEN 



- F2) * (Fl - 


F3) ) 


- F1)*(F2 - 


F3) ) 


' Fl) * (F3 - 


F2)) 


X2) .GE.O) D 


= (XI 


F, Z, BMIX) 




,) RETURN 





XI = 


X3 


Fl = 


F3 


ELSE 




X2 = 


X3 


F2 = 


F3 


ENDIF 




ELSE 





.Swap in new value from quadratic solution 
IF (F*F3.LT.0) THEN 
XI = D 
Fl = F 
X2 = X3 
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F2 = 


F3 




ELSEIF 


(F3*Fl.GT.O) 


THEN 


XI = 


D 




Fl = 


F 




ELSE 






X2 = 


D 




F2 = 


F 




END IF 






ENDIF 






CONTINUE 






D = 






END 







60 



BLOCK DATA 

INTEGER B(58) ,0(58) ,K(58) ,G(58) 

INTEGER Q{58) ,F(58) ,S{58) ,W(58) 

REAL*8 A(58) ,U(58) 

COMMON /CONSTANTS/ A, B, C, K, U, G, Q, F, S , W 

REAL* 8 MW(21) ,EI(21) ,KI(21) ,GI(21) ,QI{21) ,FI{21) ,SI{21) ,WI(21) 

REAL* 8 EIJ{21,21) ,UIJ(21,21) ,KIJ(21,21) ,GIJ(21,21) 

COMMON /PARAMETERS/ MW, EI , KI , GI , QI , Fl , SI , WI , EIJ, UIJ, KIJ, GIJ 

REAL*B XN(21) , XH(21) 
COMMON /GRENZDATA/ XN,XH 

C Equation of state parameters 

DATA A/ 
& 0.153832,600DO, 1 . 341953000D0, -2 . 998583000DO, -0 . 048312280D0 , 
& 0.375796500DO, -1 . 589575000D0, -0 . 053588470D0 , . 886594630D0 , 
& -0.710237040D0, -1 . 471722000DO , 1 . 321850350DO , -0 . 786659250D0 , 
& 0.2291290D-08, . 157672400D0 , -0 . 436386400D0 , -0 . 044081590D0 , 
& -0.003433888DO, . 032059050D0 , . 024873550D0 , . 073322790D0 , 
& -0 .001600573D0, . 642470600DO , -0 . 416260100DO , -0 . 066899570D0 , 
& 0.279179500DO, -0 . 696605100D0 , -0 . 002860589D0 , -0 . 008098836D0 , 
& 3 . 150547000DO, . 007224479D0 , -0 . 705752900D0 , . 534979200D0 , 
& -0.079314910DO, -1 . 418465000DO , -0 . 5999050D-16 , . 105840200DO , 
Sc 0.034317290D0, -0 . 007022847D0 , . 024955870D0 , . 042958180D0 , 
& 0.746545300DO, -0 . 291961300DO , 7 . 294616000DO , -9 . 936757000DO , 
& -0.005399808DO, -0 . 243256700DO , . 049870160D0 , . 003733797D0 , 
& 1.874951000D0, . 002168144D0 , -0 . 658716400DO , . 000205518D0 , 
& 0.009775195DO, -0 . 020487080D0 , . 015573220D0 , . 006862415D0 , 
& -0.001226752D0, . 002850908DO/ 

DATA B/1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,1, 1,2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 
& 3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,5,5,5,5,5,6,5,7,7,8,8,8,9,9/ 

DATA C/0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 1,1, 1,1, 1,0, 0,1, 1,1, 1,1, 1,1, 0,1, 
& 1,1,1,1,1,1,1,1,0,0,1,1,1,1,1,0,1,1,1,1,0,1,0,1,1,1,1,1,1/ 

DATA K/ 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 3, 2, 2, 2, 4, 4, 0,0, 2, 2, 2, 4, 4, 4, 4, 0,1, 
& 1,2,2,3,3,4,4,4,0,0,2,2,2,4,4,0,2,2,4,4,0,2,0,2,1,2,2,2,2/ 

DATA G/ 0,0, 0,0, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0,1, 
& 0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1,0,1,0,0/ 

DATA Q/0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 1,0, 0,0, 0,0, 0,0, 0,0, 1,0, 1,0, 
& 0,0,0,0,0,0,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,0,1,0,0,0,0,0,1/ 

DATA F/0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,1, 0,0, 
& 1,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/ 

DATA S/0, 0,0, 0,0, 0,0, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 
& 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/ 

DATA W/ 0,0, 0,0, 0,0, 0,0, 0,1, 1,1, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0, 0,0,0, 
& 0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0/ 

DATA U/ , . 5D0 , 1 , 3 . 5D0 , -0 . 5D0 , 4 . 5D0 , . 5D0 , 7 . 5D0 , 9 . 5D0 ,6,12,12. 5D0 , 
& -6,2, 3,2,2,11,-0.5D0, . 5D0 , , 4 , 6 , 21 , 23 , 22 , -1 , -0 . 5D0 , 7 , -1 , 6, 
Sc 4,1,9,-13,21,8,-0.5D0,0,2,7,9,22,23,1,9,3,8,23,1.5D0,5, 

& -0.5D0,4,7,3, 0,1,0/ 

C Characterization Parameters 

DATA MW/16 .043 0D0, 28.0135D0, 44.0100D0, 30.0700D0, 44.0970D0, 

& 18.0153D0, 34.0820D0, 2.0159D0, 28.0100D0, 31.9988D0, 

& 58.1230D0, 58.1230D0, 72.1500D0, 72.1500D0, 86.1770D0, 

& 100.2040DO, 114.2310D0, 128.2580D0, 142.2850D0, 4.0026D0, 
& 39.9480D0/ 
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DATA 


EI/151. 318300D0. 


99.737780D0, 


& 


298.118300D0, 


514.015600D0, 


& 


105 . 534800D0, 


122.766700D0,. 


& 


365.599900D0, 


370.682300D0, 


& 


450.325022D0, 


470.840891D0, 


& 


119.629900D0/ 




DATA 


KI/0.4619255D0, 


0.4479153D0, 


& 


0.5837490D0, 


0.3825868D0, 


& 


0.4533894D0, 


0.4186954D0, 


& 


0.6738577D0, 


0.6798307D0, 


& 


0.7849550D0, 


0.8152731D0, 


& 


0.4216551D0/ 




DATA 


GI/0, 0. 


027815D0, 0.18 


& 


0.332500D0, 0. 


088500D0, 0.03 


& 


0.256692DO, 0. 


,281835D0, 0.33 


& 


0.337542D0, 0. 


,383381D0, 0.42 


DATA 


QI/2*0, 0.69D0, 


2*0, 1.06775D0 


DATA 


FI/7*0, 1, 13*0/ 


1 


DATA 


SI/5*0, 1.5822D0, 0.390DG, 14* 


DATA 


WI/5*0, 1, 15*0/ 


Bina: 


ry interaction ps 


irameters 



241.960600D0, 244 . 166700D0 , 

296.355000D0, 26.957940D0, 

324.068900D0, 337 . 63 8900D0 , 

402.636293D0, 427 . 722630D0 , 

4B9.558373D0, 2.610111D0, 

.4557489D9, 0.5279209D0, 

.4618263D0, 0.3514916DO, 

.6406937D0, 0,6341423D0, 

.7175118D0, 0.7525189D0, 

.8437826D0, 0.3589888D0, 



DATA EIJ/441*1/ 
DATA UIJ/441*!/ 
DATA KIJ/441*!/ 
DATA GIJ/441*1/ 



9065D0, 0.079300DO, 
4369D0, 0.038953D0, 
2267D0, 0.366911D0, 
7354D0, 0.469659D0, 
, 0.633276D0, 14*0/ 

0/ 



0.141239D0, 
0.021000D0, 
Q.289731D0, 
0, 0/ 



DATA 



DATA 



DATA (EIJ(1, J) , J= 
0.971540D0, 
0.931484D0, 
0.989844D0, 
0.880973D0, 
(EIJ{2, J) , J= 
1.022740D0, 
1.086320D0, 
0.959340D0, 
(EIJ(3 , J) , J^ 
.925053DO, 
1.5D0, 
0.859754D0, 
0.765171D0/ 

(EIJ(4, J) , J 
1.164460D0, 

(EIJ(5, J) , J 

(EIJ(7, J) , J 
1.012821D0, 

(EIJ(8, J) , J 



^2,19) 
0.960 
1.170 



002 
881 
14) 



/ 

644D0, 

520D0, 

350D0, 

067D0, 

/ 

120D0, 

710D0, 

520D0/ 

/ 

237D0, 



0.970 
1.005 
0.945 
^4,19) 
0.960 
1, 
0.855134D0, 



1, 0.994635D0, 

0.990126D0, 1, 

0.999268D0, 1.107274D0, 
0.881161D0/ 

0.945939D0, 0.746954D0, 

1.021000D0, 0.946914D0, 



0.849408D0, 0.955052D0, 
0.906849D0, 0.897362DO, 
0.831229D0, 0.808310D0, 



0.708218D0, 
1.019530DO, 
0.880880D0, 



0.902271D0, 
0.973384D0, 



1.281790D0, 
0.726255D0, 
0.786323D0, 



DATA 

DATA 
DATA 

DATA 



:5,14) /1.022560D0, G.693168D0, 0.946871D0, 

3*1, 1.013060D0, 1, 1.00532D0/ 
:B,12) /1.034787D0, 3*1, 1.Q049D0/ 
:15,19)/1.008692D0, 1.010126D0, 1.011501D0, 

1.014089D0/ 
:9,12) /I.IDO, 1, 1.3D0, 1.3D0/ 



DATA (UIJ(1, J) , J=2, 

& 0.886106D0, 0. 

& 1.156390D0, 3' 

& 1.191904D0, 1. 

DATA (UIJ(2, J) , J=3, 

& 0.993476D0, 0. 

DATA {UIJ(3 , J) , J=4, 



19) / 

963827D0, 1, 0.990877D0, 1, 0.736833D0, 
1, 0.992291D0, 1, 1.003670D0, 1.302576D0, 
205769D0, 1.219634D0, 1.233498D0/ 



12) /0.835058D0 
408838D0, 3*1, 
19) /0.969870D0 
1.077634D0, 



& 5*1, 1.066638D0 
& 1.108021D0/ 

DATA (UIJ{4, J) , J=5,14) / 
& 1.065173D0, 1, 0.971926D0 

DATA (UIJ(7, J) , J=7,19) /8*1, 
Sc 1.038338D0, 1.042735D0 



0.816431DO, 0.915502D0, 1, 
0.993556D0/ 

, 2*1, 1.045290D0, 1, . 9D0 , 
1.088178D0, 1.098291DO, 



.616560D0, 2*1, 4*1.25D0/ 
1.028973D0, 1.033754DO, 
1.046966D0/ 



DATA (KIJ{1, J) , J=2, 19) / 
& 1.003630D0, 0.995933DO, 1, 1.007619D0, 1, 1.000080D0, 
& 1.023260D0, 3*1, 0.997596D0, 1, 1.002529D0, 0.982962D0, 
& 0.983565D0, 0.982707D0, O.981849D0, 0.980991D0/ 

DATA (KIJ(2 , J) , J=3 , 8) / 
& 0.982361D0, 1.007950D0, 1, 1, . 9425960D0 , 1 . 032270DO/ 

DATA (KIJ(3,J) , J=4,19)/1.008510D0, 2*1, 1 . 00779D0 , 7*1 . ODO , 
& 0.910183DO, 0.895362D0, 0.881152D0, 0.867520D0, 0.854406D0/ 
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DATA (KIJ(4, J) , J=5,8) /0.986893D0, 1, 0.999969D0, 1.020340D0/ 
DATA (KIJ(7, J) , J=7,21) /8*1, 0.958130D0, 0.962870D0, 
& 0.957828D0, 0.952441D0, 0.948338D0, 2*1/ 

DATAGIJ{1,3) /0.807653D0/ 

DATAGIJ(1,8) /1.957310D0/ 

DATAGIJ{2,3) /0.982746D0/ 

DATAGIJ(3,4) /0.370296D0/ 

DATAGIJ(3,6) /1.673090D0/ 



DATA XN/ 






















& 50.0 , 


0.0 , 


0.0 , 


0.0 , 


0.0 , 


0, 


.0 , 


0, 


.0 , 


0, 


, 


& 0.0 , 


0.0 , 


0.0 , 


0.0 , 


. 0.0 , 


0, 


.0 , 


0, 


.0 , 


0, 


. 


& 0.0 


0.0 , 


0.0 , 


0.0 , 


0.0/ 














DATA XH/ 






















& 100.0 , 


20.0 , 


5.0 , 


1.5 , 


, 1.5 , 


0, 


.5 , 


. 0, 


.5 , 


0, 


.1 


& 0.05 , 


0.05 . 


0.05 , 


0.05 , 


, 30.0 , 


50, 


.0 , 


, 0. 


.02 , 


. 


.5 


& 0.015, 


0.02 , 


0.02 , 


10.0 , 


, 3.0/ 














END 























27 



IS 15305 {Part 2) : 2003 
ISO 12213-2 : 1997 

Annex G 

(informative) 

Bibliography 

[1] Starling, K.E., Savidge, J.L. "Compressibility Factors for Natural Gas and Other Related Hydrocarbon 
Gases", American Gas Association (AGA) Transmission Measurement Committee Report No. 8, American 
Petroleum Institute (API) MPMS, chapter 14.2, second edition, November 1992. 

[2] Jaeschke, M., Humphreys, A.E. "The GERG Databank of High Accuracy Compressibility Factor 
Measurements", GERG Technical Monograph TM4 (1990) and Fortschritt-Berichte VDI, Series 6, No. 251 
(1990). 

[3] Jaeschke, M., Hinze, H.M., Humphreys, A.E. "Supplement to the GERG Databank of High Accuracy 
Compressibility Factor Measurements", GERG Technical Monograph TM7 (1996) and Fortschritt-Berichte VDI, 
Series 6, No. 355 (1996). 

[4] SCHOUTEN, J.A., MiCHELS, J.P.J. "Evaluation of PUT Reference Data on Natural Gas Mixtures — Final report", 
Appendix to Gas Research Institute Report No. GRI/93-006, September 1992. 

[5] Savidge, J.L., BEYERLEiNf, S., Lemmon, E. Technical reference document on the 2nd edition of AGA Report 
No. 8, November 1992 (Gas Research Institute Report No. GRi/93-0181, May 1993). 



28 



Bureau of Indian Standards 

BIS is a statutory institution established under the Bureau of Indian Standards Act, 1986 to promote 
harmonious development of the activities of standardization, marking and quality certification of 
goods and attending to connected matters in the country. 

Copyright 

BIS has the copyright of all its publications. No part of these publications may be reproduced in any 
form without the prior permission in writing of BIS. This does not preclude the free use, in the course 
of implementing the standard, of necessary details, such as symbols and sizes, type or grade 
designations. Enquiries relating to copyright be addressed to the Director (Publication), BIS. 

Review of Indian Standards 

Amendments are issued to standards as the need arises on the basis of comments. Standards are also 
reviewed periodically; a standard along with amendments is reaffirmed when such review indicates that 
no changes are needed; if the review indicates that changes are needed, it is taken up for revision. 
Users of Indian Standards should ascertain that they are in possession of the latest amendments or 
edition by referring to the latest issue of 'BIS Catalogue' and 'Standards: Monthly Additions'. 

This Indian Standard has been developed from Doc: No. PCD 24 (2031). 

Amendments Issued Since Publication 



Amend No. 



Date of Issue 



Text Affected 



BUREAU OF INDIAN STANDARDS 
Headquarters: 

Manak Bhavan, 9 Bahadur Shah Zafar Marg, New Delhi 1 10002 
Telephones: 2323 0131 , 2323 3375, 2323 9402 

Regional Offices: 

Central : Manak Bhavan, 9 Bahadur Shah Zafar Marg 
NEW DELHI 110002 

Eastern : 1/14 C.I.T. Scheme VII M, V.I.R Road, Kankurgachi 
KOLKATA 700054 

Northern : SCO 335-336, Sector 34-A, CHANDIGARH 160022 
Southern : C.I.T. Campus, IV Cross Road, CHENNAI 6001 13 



Western 



Branches 



Manakalaya, E9 MIDC, Marol, Andheri (East) 
MUMBAI 400093 



Telegrams: Manaksanstha 
(Common to all offices) 

Telephone 

(-2323 7617 
12323 3841 

r2337 8499, 2337 8561 
12337 8626,2337 9120 

(-60 3843 
160 9285 

r2254 1216, 2254 1442 
12254 2519,2254 2315 

r 2832 9295, 2832 7858 
I 2832 7891 , 2832 7892 



AHMEDABAD. BANGALORE. BHOPAL. BHUBANESHWAR. COIMBATORE. FARIDABAD 
GHAZIABAD. GUWAHATI. HYDERABAD. JAIPUR. KANPUR. LUCKNOW. NAGPUR. 
NALAGARH. PATNA. PUNE. RAJKOT THIRUVANANTHAPURAM. VISAKHAPATNAM. 



Printed at Simco Printing Press, Delhi 



